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ABSTRACT 

Aims. We present a quantitative study of a new data set of high redshift Type la supernovae spectra, observed at the Gemini telescopes 
during the first 34 months of the Supernova Legacy Survey. During this time 123 supernovae candidates were observed, of which 87 
have been identified as SNe la at a median redshift of z = 0.720. Spectra from the entire second year of the survey and part of the 
third year (59 total SNe candidates with 46 confirmed SNe la) are published here for the first time. The spectroscopic measurements 
made on this data set are used determine if these distant SNe comprise a population similar to those observed locally. 
Methods. Rest-frame equivalent width and ejection velocity measurements are made on four spectroscopic features. Corresponding 
measurements are presented for a set of 167 spectra from 24 low-z SNe la from the literature. 

Results. We show that there exists a sample at high redshift with properties similar to nearby SNe. The high-z measurements are 
consistent with the range of measurements at low-z and no significant difference was found between the distributions of measurements 
at low and high redsift for three of the features. The fourth feature displays a possible difference that should be investigated further. 
Correlations between Type la SNe properties and host galaxy morphology were also found to be similar at low and high z, and within 
each host galaxy class we see no evidence for redshift-evolution in SN properties. A new correlation between SNe la peak magnitude 
and the equivalent width of Sill absorption is presented. Tests on a sub-set of the SNLS SNe demonstrates that this correlation reduces 
the scatter in SNe la luminosity distances in a manner consistent with the lightcurve shape-luminosity corrections that are used for 
Type la SNe cosmology. 

Conclusions. We show that this new sample of SNLS SNe la has spectroscopic properties similar to nearby objects. 
Key words, supernovae: general - cosmology: observations - surveys 

1. Introduction |2004|) . but the results are still consistent with a wide range 

, . n , . of dark energy models. Th e CFHT Supernova Legacy Survey 

Hubble diagrams m akmg use of die standardizeable na - rSNLS.' Astier all (l2006h l aims to increase the measured pre- 



tore of Typ e laS Ne (|Riess a . 11 19981 ; | Perlmuttere?a/. || l999t ^-^-^^ ^^-^ parameter by finding and following about 500 dis 

\MmeLaimm Astwr et alA |20^ indicate the significant ^^^^ ^ ^^^^ ^ ^ ^-^^^-^^ 

presence of a dark energy that is driving the accelerating ex- 

pansion of the universe. This result is additionally constrained The fundamental assumptions of the methods used to es- 

by observati ons of the pow er spectrum of the cosmic microwave Annate the cosmological parameters with SNe are that the se- 

backrou nd (ISpergel e/a/.l [2 003. 2006), properties of massive l^'^^ed low-z and high-z objects have similar peak magnitodes 



clusters ('Alie^n.Schimdt & Fabian 2002) and other large scale ^1*^" l^'^^H^ observed dispersion) and that the es- 
structure observations ( Seliak et al. 2005). The equation of state tablished hghtcurve shape-luimnosity relations are applicable 
parameter of this dark energy (defined as w = /,/p, the ratio of ^° *ese SNe at all redshifts. However, the current under- 
dark energy pressureptodensitvp)haste^ standing of SNe la explosions [see Hdlebrandt & NiemeyeiJ 
recent surveys (iKnop et al\\200± (Barris a/.lllool lRiess7fa/.l <MD a review] and their progenitors (PGs) illustrates that 



there may be some differences, or evolution, in their inherent 
Send offprint requests to: 1. J. Bronde^r properties at high redshifts (Hoflich et al 1998, Hatano et al 

* The Figures in Appendix A are only available in electronic form via 2000, Lentz et al 2000). Such effects could distort the cos- 
|http://www.edpsciences/anada.org mological results of SNe la surveys. The differences in prop 
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erties of SNe in hosts o f different morpholo gy |Wang et aU 
119971 lUmedaef [1991 iHamuv et al.\ UWOl lHowelL " ,200lt 
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ISuUivan et a/1l2003l: iGallagher et alll2005t ISuUivan et a/ll2006l) 
provides additional evidence of SNe la heterogeneity. The 
luminosity-lightcurve shape calibrations used for SNe l a cos- 
mology (Riess, Press, & Kirshner 1996; iPhillips et all 119991; 
iGoldhaber et aZ..2001;,Guv et aL200 5) make an empirical acco- 
modation for the dispersion of SNe properties observed locally, 
but these methods may be vulnerable to systematic changes in 
the SNe la population at high redshifts (Leibundgut 2001 ). Thus, 
it is imperative to undertake a quantitative investigation to deter- 
mine whether there are any changes in these objects that may 
affect their reliability as distance indicators for cosmology. 

The spectra of SNe la reveal more subtle characteristics than 
photometric observations and are well suited for this investiga- 
tion. For example, measurements of photospheric ejection ve- 
locities (vej) explore the kinetic energy and distribution of el- 
ements in the expanding SN ejecta. These measurements pro- 
vided some o f the first evid ence of intrinsic differences between 
Type la SNe dBranch et al. il988,) and have also been shown t o 



correlate with luminositv (l\yells et al. ] ll994HFisher ef a/.lll995h . 



methods for measureing EW and Vej are presented in Sect. [3] 
The results and a comparison of the two SNe sets are presented 
in Sect. |4l in Sect. |5] these results are discussed in more detail 
and the conclusions are presented in Sect. |6l 



2. Spectroscopic Data Sets — Overview 

A large number of the SNe la that were used to provide 
the latest constrain t s on w and other co smological parameters 
dAstier et al. Il2006t ISpergel et a/.ll2006l) are studied in this pa- 
per. Specifically, this includes the SNLS objects observed at the 
Gemini telescopes during the first year of the survey. This set is 
augmented by Gemini spectra from the entire 2"'' year and part 
of the year of SNLS observations, which are published here 
for the first time. This section briefly describes the observation, 
data reduction, and identification of these SNLS SNe spectra in 
Sect. 12.11 In Sect. 12.21 a summary of the set of low-z SNe that 
were quantitatively compared to the distant SNe la is provided. 



More recentlv, [Benetti et al\ i. 2005 ^ used cluster analysis of mea 

surements of the decline rate of SNe la Vej to yield a new way of ^.1. High Redshift SNLS Supemovae 



grouping low-z SNe that perhaps illuminates the role of differ- 
ent explosion mechanisms in SNe la diversity. Comparisons of 
SNe la ejection veloci ties from ob j ects a t different redshifts hav e 
also been made [e^. iHook et al.\ (l2005h ; iBlondin et al\ (l2006h ; 
^allandeffl/. (20061)] demonstrating that the range of veloci- 
ties measured in SNe la ejecta are comparable at all redshifts. 
Folatelli (2004, F04 hereafter) illustrated the utility of equiv- 
alent width (EW) me asurements on low-z SNe la spectra and 
iGaravini et al.\ (l2007h extended this work to compare 12 high-z 
object s to the low-z population. H achin ger, Mazzali, & Be nettii 
(120061) also illustrated correlations between these EW mea- 
surements an d the vel ocity-gradient clas sification scheme from 
iBenetti et a/.l (2005). B ranch et al.\ ( l2006h presents additional ev- 
idence for a temperature or explosion mechanism sequence in 
SNe la with comparisons of different low-z objects in EW-space. 

These studies all illustrate how quantitative spectroscopic 
measurements can be used to explore the factors that affect the 
heterogeneity and homogeneity of SNe la. The utility of these 
parameters for exploring differences in high-z SNe la is limited, 
however, as the majority of these measurements cannot be repro- 
duced in the high-z data. These limitations are due to the reduced 
signal-to-noise (S/N), redshift effects (which shift the features to 
redder regions where sky noise is dominant or moves features 
out of optical detection limits entirely), host galaxy contamina- 
tion and the smaller number of observations which all charac- 
terise high-z data. 

The analysis presented here includes rest frame equivalent 
width and ejection velocity measurements on a large set of SNe 
la spectra observed for the SNLS with methods that are specif- 
ically tailored for high-z data. Th ese measurements are similar 
to the recently published work of IGaravini et al.\ (l2007h but are 
utilized on a much larger high-z set and include a different treat- 
ment of host galaxy contamination and estimates of the low-z 
trends. The studied set of SNe were selected from the objects 
considered 'confirmed' SNe la by the SNLS [see Howell et al. 
(2005), H05 hereafter] and used in the survey's cosmological 
analysis. These high-z results are compared to corresponding 
measurements on a set of low-z supemovae from the literature. 
The results of this investigation address the cosmological impli- 
cations from high-z SNe surveys by determining whether this 
SNLS-set of distant objects comprises a population that is spec- 
troscopically similar to their low-z counterparts. An overview of 
the data sets used in this analysis is presented in Sect. |2l The 



2.1 .1 . SNLS Spectra - Observation, Data Reduction and 
Classification 

Full details of the methods used for SN LS SN candidate se- 
lection and observation can be found in Sullivan et al. (20061) 
and H05, respectively. In summary, the Canada-France-Hawaii 
Telescope Legacy Survey (CFHTLS), using the Megacam wide 
field imager, acquires real time lightcurves of possible SNe by 
repeatedly imaging four, one square-degree fields approximately 
every four days (observers' frame) in several filters (a combina- 
tion of g', r', i', and z') as part of the Deep component of the 
CFHTLS. Spectroscopic observations are then made at 8m-class 
telescopes in order to identify the SN-type and redshift of these 
candidates. 

The objects studied here were observed at the Gemini 
North and South telesc opes with the G emini Multi Object 
Spectrograph [GMOS; dHook et al.\ l2004 l in longslit mode. 
Typically, three 20 minute to 30 minute exposures are taken of 
each candidate. The GMOS R400 grating (400 lines per mm) 
and 0.75" slit are used, which gives a 6.9 A spectral resolution. 
A central wavelength of 680nm was chosen for objects at an es- 
timated z < 0.4 and objects at higher estimated z are observed 
with a central wavelength setting of 720nm. This central wave- 
length selection is important in order to ensure that the defin- 
ing Sill features (located at rest wavelengths of ~ 4100 A and 
~ 6150 A) of SNe la spectra are observed. The 680 nm setting 
was paired with the GG455_G0305 order-sorting filter to provide 
wavelength coverage from 4650 A to 8900 A. The 720 nm set- 
ting was used along with the OG515_G0306 order-sorting filter 
to cover wavelengths from 5100 A to 9300 A. Queue observing 
was employed to observe the SNLS candidates, which require 
better than 0.75" image quality (corrected to zenith), in either 
classical or nod-and-shuffle (N&S) observing modes (see H05). 

The reduction of the spectra was completed using 
independently from that of H05, but with very similar methods. 
The independent reductions were made to check for any system- 



' IRAF is written and supported by the IRAF programming group 
at the National Optical Astronomy Observatories (NOAO). NOAO is 
operated by the Association of Universities for Research in Astronomy 
(AURA), Inc. under cooperative agreement with the National Science 
Foundation 
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atic differences in the pipelines and the results from all reduc- 
tions were in very close agreement. 

To classify the SNLS objects, the final spectra are quanti- 
tatively compared to template supernovae fro m large library of 
low-z objects with ax^ matching program [see lHowell & Wana 
(l2002h and H05 for details]. The results fr om this matching pro- 
gram and the opinions of SNLS members are used to rank the 
identifications according to one 'confidence index' (CI) from the 
six categories defined in H05, with a 5 being a 'Certain la' and 
a as 'Not a la' . Objects that are ranked at a CI of 3 or greater 
are considered 'confirmed' Type la SNe. These are the objects 
that were used in the cosmology fits from lAstier et al. I (l2006l) 
and also make up the 80 high-z SNe la studied here. 

The selection criteria of the SNLS sample means that Type la 
SNe at high-z with properties that are very different from those 
seen nearby are likely to be rejected. Thus in the comparisons 
below the existence of a new population of high-z SNe with un- 
usual properties cannot be ruled out. However as will be seen 
later, the selection process does allow a range of SNIa proper- 
ties through, and thus tests for evolution within this range are 
possible. 

2.1.2. Summary of the SNLS Dataset 

The SNe candidates observed during the first year of the survey 
at the Gemini telescopes were published in H05. Table 2 of H05 
lists the derived properties of these spectra. This set includes 64 
candidates observed from August 2003 through October 2004. 
Within this set, 41 objects were classified with as confirmed SNe 
with a CI > 3 and are studied here. 

Spectra from the entire second year and part of the third year 
of SNLS observations at Gemini are presented here for the first 
time. These spectra were observed between November 2004 and 
May 2006. During this period 59 SNe candidates were observed 
and 46 of these objects were confirmed as SNe la. The setup and 
conditions of the observations of these objects are listed in Table 
|6l The derived properties and final classifications of these spectra 
are summarized in Table|7] Some example spectra are presented 
in Fig.[T]and all of the spectra are available in the on-line version 
of this paper 

For the analysis in this paper we make use of not only spec- 
tral properties but also photometric properties. The latter are de- 
rived from lightcurve fits to the SNLS CFHT photo metry, fol- 
lowing the methods described in lAstier et al. I (l2006l) . The main 
lightcurve parameters used are the date of B-band maximum 
light, and percentage increase relative to the host galaxy (inter- 
polated to the date of spectroscopy). 



2.2. Low Redshift Type la Supemovae 
2.2.1. Low-z Spectra 

The spectra that make up the low redshift sample in this study 
were selected from the literature. A majority of these spec- 
tra can be found on the online Supernova Spectrum Archive 
(SUSPECT0)- The main goal of this study is to quantify the sim- 
ilarities and differences between high-z SNe la and those dis- 
covered locally, so objects with epoch and wavelength coverage 
that approximated that of the SNLS sample were selected. These 
spectra were corrected for peculiar velocities and redshift where 
appropriate. All of the low-z objects that met this criteria are 
listed in Table [T] 



The low-z sample includes 167 spectra from 2 4 different ob- 
jects. The set includes the so-called 'core normal' (Branc h et al. 
2006) objects [i.e. SNe si milar to SN 1981B (Branch ef a, 
Il983h . SN 1989B (Barbon ^7ZlfT990 : Wells ef a/.lll994l) . SN 



http : //bruf ord . nhn . ou . edu/suspect/ index 1 . html 



1992A jKirshner et a/...1993i ror SN 1994D (.Patat ef q/.l 1 19961) 1 
as well as spectroscopically peculiar SNe la such as the overlu- 
minous 1 99 IT-like objects and underluminous 1991bg-like SNe. 

2.2.2. Low-z photometry and distance estimation 

The photometric properties of these nearby objects are listed 
in Table |2] T hese propertie s have been presented in the litera- 
ture (.Perlmutter ef 'a/.lfl999l: iPhillips et al.\\l99§i. lAltavilla et al.\ 
120041: iReindl et a/"1l2005h with slightly different results depend- 
ing on the color corrections used or the lightcurve fits to the 
published photometry. To remove systematic differences when 
comparing results, the photometry from these SNe was re-fit 
by SNLS members to derive the magnitudes and lightcurve 
shape parameters used later in this analysis. This lightcurve 
fitti ng used a new te mplate spectrum for K-corrections [see 
iHsiao et al. (in press^ and Conley et al. (in prep) for full details]. 
The parameters derived from these fits are summarized in Table 

m 

These re-derived photometric properties were also used to 
estimate the absolute peak (B') magnitudes, M^^^^^, of the low- 
z SNe la. A new estimate of these magnit udes is of interest as 
the published values hav e [in general, see iPhillips et al.\ (1 19991) 
andl Altavilla et a/.l(l2004 l) for exceptions] already been corrected 
for a derived lightcurve shape-luminosity correction. In order to 
fully explore the homogeneity of SNe at different redshifts and 
to investige any possible spectroscopic sequences in SNe la, an 
uncorrected, or intrinsic, M^^,^ is desired. 

To make this estimation, the observed magnitudes from 
Table |2] were first corrected for host galaxy reddening using the 
relation between Amis(B) (estimated from the re-fit ligh t curves ) 
and B - V colour described by Eq. 7 in iPhillips et al.\ ( Il999l) . 
Combined with the relations in Eq. 9 and Eq. 1 1 from the same 
publications, the observed and estimated 'true' B - V colours 
were approximated to constrain the host galaxy reddening for 
each SN. This is the only step of the estimate that made use of 
a SNe la lightcurve shape-luminosity correction and is unavoid- 
able as some assumption of the intrinsic colour of these objects 
is necessary to infer the host galaxy reddening. This reddening 
was then converted to a B' -extinction with Rb = 3.0. This is not 
the standard value used in this exercise and does not take into 
account extremely red objects ( Elias-Rosa et a/.ir2006 ^. but it is 
approximately the median value of the different Rb estimates in 
the literature. The selection of this Rb value does not have a large 
effect on the results in this paper. A different Rb would systemat- 
ically shift all the values in Table|2]a slight amount, but leave the 
relative relationship between the magnitude of the different SNe 
the same, which is the true parameter of interest in this study as 
illustrated in Fig.fTOl 

Next, the distance modulus, jj., to the host galaxy was used 
to convert to an absolute magnitude via fi - m - M (where 
m and M are the observed and absolute magnitudes, respec- 
tively). To yield the most correct distance modulus, the authors 
used the most reliable Cepheid variable based distances to the 
SNe host galaxies that were available in the literature. These 
distances were available for 9 of the SNe in Table [T] and an- 
other 5 low-z SNe had reliable host distances from the Nearby 
Galaxies Catalogue (Tullv 1988). The remaining objects did not 
have any independent distance constraints and were not consid- 
ered for these Mb estimates. The errors in each of these steps 
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Fig. 1. Example spectra of SNLS candidates observed by Gemini. The observed spectra after subtraction of the best-fitting host 
galaxy template are displayed in green. The smoothed spectra are shown in red and the best fitting template SNe are overplot in 
black. The top two objects are 'certain type la' SNe (CI - 5) as they have clear signs of Sill. The bottom objects are SNLS candidates 
with CI = 4 (left) and CI = 3 (right). These spectra match SNe la templates but have ambiguous Sill absorption. All of the other 
spectra are available online. 



Table 3. Wavlength regions for the SNela features defined by 
F04. 



Feature 


Label 


Blue-ward 


Red-ward 


Number 




Limit (A) 


Limit (A) 


1 


'Call H&K' 


3500 - 3800 


3900 - 4100 


2 


'Sill 4000' 


3900 - 4000 


4000 - 4150 


3 


'Mgir 


3900 - 4150 


4450 - 4700 


4 


'Feir 


4500 - 4700 


5050 - 5550 


5 


'SII W 


5150-5300 


5500 - 5700 


6 


'Sill 5800' 


5550 - 5700 


5800 - 6000 


7 


'Sill 6150' 


5800 - 6000 


6200 - 6600 



were propagated through this estimate to provide a final error on 
Mb. The final column in Table|2]lists these intrinsic peak magni- 
tude values for the 14 low-z SNe that had reliable, independent 
distance estimates and secure observed photometric properties. 

3. Spectroscopic Analysis Techniques 

3.1. Equivalent Widths 

Rest frame equivalent width measurements are a shape- 
independent method of quantifying the strength of spectroscopic 
feat ures. The method use d in this EW analysis is similar to F04 
and lGaravini et al. I (l2007h . which defined eight regions that con- 
tain the strongest absorption features observed in core normal 




0.0 r . I . . . . I , , , . I . . . . I , , , , I . . . . I , , , 

3500 4000 4500 5000 5500 6000 6500 
rest wovelength (Angstroms) 



Fig. 2. The spectroscopic regions of specific interest to SNe la 
EW measurements, as defined by F04, are shown here on a core 
normal SN la near maximum light. The wavelength ranges and 
nomenclature for these features is summarized in Table [3] 



SNe la spectra near maximum luminosity. The first seven regions 
are illustrated in Fig. |2] and the nomenclature and wavelength 
ranges used to define them are given in Table [3] For supernovae, 
the continuum needed to make the EW measurement is not well 
defined, so it is estimated by making a straight line fit between 
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Table 1. The low redshift spectra analyzed in this study. The dates correspond to the epoch relative to 
maximum luminosity in the B-band. Many of these spectra are available in the SUSPECT archive as well 
as their listed publications. 



Name 


Date(s) 


Reference 


1981B 


0,20 


Branch et al. 1983 


1986G 


-5,-4,-3,0 


Phillips et al. 1987 


1989B 


-7,-5,-3,-2,-1,3,5,8,9,11,12,13,16,17,18,19 


Barbon et al. 1990, Wells et al. 1994 


1990N 


-14,-7,2,3,17 


Mazzali et al. 1993 


1991M 


3 


Gomez & Lopez 1998 


1991T 


-12,-11,-10,-8,-7,-6,0,10,15 


Mazzali, Danziger, & Turatto 1995 


1991bg 


0,1,2 


Turattoe/a/. 1996 


1992A 


-5,-1,3,5,6,7,9,11,16,17 


Kirshner ef a/. 1993 


1994D 


- 1 1,-10,-9,-8,-7,-5,-4,-2,2,4, 10, 1 1 , 1 2, 1 3,20 


Patat et al. 1996 


1994Q 


10 


Gomez & Lopez 1998 


1996X 


-2,-1,0,1,2,7 


Salvo et al. 2001 


1997br 


-9,-8,-7,-6,-4,8 


Li et al. 1999 


1997cn 


3 


Turatto e/ a/. 1998 


1998aq 


-8,0,1,2,3,4,5,6,7 


Branch et al. 2003 


1998bu 


-4,-3,-2,8,9,10,11,12,13 


Capellaro et al. 2001 


1999aa 


-11,-3,-1,5,6 


Garavini et al. 2004 


1999ac 


-15,-9,0,2,8,11 


Garavini et al. 2004 


1999aw 


3,5,5,12 


Strolger et al. 2002 


1999by 


-5,-4,-3,-2,3,4,5,6,7,8,11 


Garnavich et al. 2004 


1999ee 


-1 1,-9,-8,-6,-4,-2,- 1 ,0,3,5,7,8,9, 12, 14 


Hamuy et al. 2002 


2000E 


-2,3,5,14 


Valentini et al. 2003 


2000cx 


-4,-3,-2,-1,0,1,5,6,7,9,11 


Li etal. 2001 


2002bo 


-14,-13,-11,-7,-6,-5,-4,-3,-2,-1,4 


Benetti et al. 2004 


2003du 


-11,-7,13 


Anupama, Sahu, & Jose 2005 



the local maxima that bound each feature. The following rules 
govern how the continuum is set: 

- the pseudo-continuum is always set within the feature re- 
gions listed in Table [3] 

- the feature boundaries are selected to maximise the wave- 
length span of the measurement without intersecting the 
boundaries of any neighboring features 

- the EW is always calculated in the rest frame 

The continuum is set automatically but a human judgement 
is made on its placement to ensure the guidelines above are fol- 
lowed. The EW is then calculated by numerically integrating un- 
derneath this 'pseudo-continuum' at every wavelength pixel Aj 
(for all points within the set continuum bounds): 

where is the flux level of the spectrum, and fc is the flux of the 
pseudo-continuum. The statistical uncertainty in the EW mea- 
surement is calculated by propagating the estimated uncertain- 
ties in the flux and pseudo-continuum: 

here erf is the measurement uncertainty in the flux and ctc^ is the 
uncertainty in the pseudo-continuum flux. The 0-^. is calculated 
via the propagation of the correlated errors in the selection and 
calculation of this continuum. 

In addition to the error calculated in Eq. |2] uncertain- 
ties can be caused by variations in the extent or the slope of 
the pseudo-continuum. The range and effects of these possible 



pseudo-continuum variations were estimated by randomly shift- 
ing the continuum boundaries within | A of the selected max- 
ima (where R is the range of the feature boundaries, see Fig. 
|2]l for the feature and spectrum in question. The EW was then 
re-measured, and the standard deviation of these measurements 
(which were repeated 50 times) was accepted as an approxima- 
tion of this error in each EW measurement. For the high S/N 
spectra, the possible variations in the continuum were found to 
be the dominant source of error The high-z data was subject to 
some additional systematic errors, so this general procedure had 
to be amended for these objects; this is discussed in the follow- 
ing section. 

In calculating the low-z EW, any possible effects from noise 
fluctuations were mitigated by setting the continuum end-points 
to an average flux value that was computed over a 25 A win- 
dow around the selected maxima. Observed error spectra were 
not available for most of the objects in this nearby sample, so 
the statistical error could not be calculated. However, the dom- 
inant error from the possible pseudo-continuum variations was 
estimated as described above. In total, these uncertainties were 
on the order of 5 - 10% of the measured low-z EW. This mea- 
surement tec hniqe differs shghtly from previously published EW 
studies (FO?;'Hachinger. M azzah. & Benettill2006t [Branch et a/1 
I2OO6), but the results here are in very good agreement with these 
other references. Based on the wavelength coverage available in 
the high-z spectra, this EW analysis is focused on EW features 
1, 2, and 3 - Call H&K, Sill 4000, and Mgll in the F04 nomen- 
clature (EW{CaII), EW{SiII), & EW{MgII} hereafter). 

3.2. High-z EW Measurements 

As explained in Sect. 13.2.11 the high-z spectra were median fil- 
tered in order to set the feature bounds. The EW measurement 
was then made with the initial, un-filtered points as Fig.|3]shows. 
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Table 2. The photometric properties of the low-z spectra analyzed in this study. The published photometry 
from these objects was re-fit to determine the magnitudes and lightcurve shape parameters shown here. SN 
199 IM did not have published photometry in B-band, so the parameters for this object were not determined. 
The SNLS models did no t fit the sub-luminous SN 1991bg with any accuracy; the values given here for this 
unusual object are from Altavilla et al. \ (12004 . The precision of the derived values for the unusual SN 
2000cx are also in doubt and were not used in any of the analysis in this paper. Uncertainties in each 
quantity are given in brackets. 



Name 






ml 


(B - V)^,„ 






Mr 




1981B 


0.0060 


44671.84(0.09) 


11.83(0.01) 


-0.06(0.01) 


0.95(0.01) 


-19, 


,58(0, 


,22)*" 


1986G 


0.0027 


46560.53(0.05) 


11.91(0.01) 


0.83(0.01) 


0.73(0.01) 








1989B 


0.0024 


47564.39(0.27) 


12.12(0.02) 


0.30(0.02) 


0.94(0.01) 


-19, 


,03(0, 


,31)<" 


1990N 


0.0034 


48081.92(0.03) 


12.53(0.01) 


-0.09(0.01) 


1.07(0.01) 


-19, 


,59(0, 


23)<'> 


1991M 


0.0072 
















1991T 


0.0060 


48374.25(0.04) 


11.30(0.01) 


-0.04(0.01) 


1.05(0.01) 


-20, 


,07(0, 




1991bg 


0.0030 


48604.50(1.0) 


14.75(0.04) 


0.00(0.08) 


.68(0.10) 








1992 A 


0.0063 


48639.71(0.02) 


12.41(0.01) 


-0.09(0.01) 


0.85(0.01) 


-18, 


,73(0, 


,22)<'> 


1994D 


0.0015 


49431.41(0.01) 


11.64(0.01) 


-0.19(0.01) 


0.84(0.01) 


-19, 


,44(0, 


,22)<" 


1994Q 


0.0294 


49494.35(0.90) 


16.16(0.09) 


-0.14(0.04 


1.11(0.04) 








1996X 


0.0070 


50190.65(0.10) 


12.86(0.01) 


-0.12(0.01) 


0.90(0.01) 


-19, 


,38(0, 


,22)<" 


1997br 


0.0053 


50558.57(0.03) 


13.23(0.02) 


0.12(0.01) 


0.89(0.01) 








1997cn 


0.0162 


50593.36(0.32) 


18.12(0.02) 


0.92(0.01) 


0.89(0.02) 








1998aq 


0.0037 


50930.16(0.02) 


12.18(0.01) 


-0.25(0.01) 


0.98(0.01) 


-19 


,87(0, 


,23)<*' 


1998bu 


0.0030 


50951.89(0.04) 


11.97(0.01) 


0.15(0.01) 


0.99(0.01) 


-19, 


,18(0, 


,26)<'> 


1999aa 


0.0144 


51231.82(0.03) 


14.59(0.01) 


-0.18(0.01) 


1.12(0.01) 








1999ac 


0.0095 


51250.07(0.08) 


13.99(0.01) 


-0.08(0.01) 


1.02(0.01) 








1999aw 


0.0380 


51253.34(0.08) 


16.59(0.01) 


-0.17(0.01) 


1.22(0.01) 








1999ee 


0.0114 


51468.89(0.01) 


14.69(0.01) 


0.11(0.01) 


1.05(0.01) 


-19, 


,70(0, 


,30)<3' 


2000e 


0.0047 


51576.44(0.03) 


12.68 (0.01) 


0.01(0.01) 


1.03(0.01) 


-20, 


,08(0, 


,25)<" 


2000CX 


0.0079 


51752.03(0.03) 


12.92(0.01) 


-0.09(0.01) 


0.91(0.01) 








2002bo 


0.0042 


52356.43(0.04) 


13.82(0.01) 


0.29(0.01) 


0.95(0.01) 


-19, 


,04(0, 


,46)<^> 


2003du 


0.0064 


52766.22(0.19) 


13.36(0.01) 


-0.25(0.01) 


1.02(0.01) 


-19, 


,53(0, 


,45)<5> 



a ~ date of maximum luminosity (rest-frame B') 

b - maximum apparent magnitude in rest-frame B' magnitudes, corrected for Milky Way extinction with E(B-V) values 

from Schlegel et al. (1998) 

c - 5 - V color at (B-band) maximum luminosity 

d - stretch values, based on fitting a fiducial s = 1.0 Type la template 

e — absolute peak luminosities (B'), not adjusted for lightcurve shape-luminosity corre ctions. These made use of 
the following distance modulus references - 1 : Gibson efaZ. (2000) 2: Sahae/a/. ( 20 0i|) 3: IStritzinger et al\ ( |2002|) 4: 
ICandia et fl/.l l l2003b 5: lAnuDama. Sahu. & Josel ( 120051) 6: IRiess et a/.l ( 120051) 7: lTullvi ( ll^l) 

3.2.1 . High-z EW Systematic Errors I - Signal to Noise 



The sky noise affects the high-z measurements by augmenting 
the local maxima selected to set the pseudo-continuum. This can 
result in an overestimate of the equivalent width. A median filter 
was used to correct this systematic effect. The size of the filter 
window was selected via simulations in which spectra from the 
low-z sample were degraded with noise to match that in the high- 
z SNe. The EW was then re-calculated on each of the artificial 
spectra, using a median filter set to various window sizes. Using 
a median filter over a 60 A window proved to be the most robust 
for all measurements at the wide range of S/N levels covered by 
the high-z spectra. 



These simulations were also used to check that the uncer- 
tainties in EW measurements behaved as expected for a given 
signal-to-noise input spectrum. Both the simulations and the fi- 
nal measurement results demonstrated that the statistical errors 
in the high-z sample were approximately 7% - 20% of the mea- 
sured EW value. This is nearly twice the uncertainty of the low-z 
objects due to their reduced S/N. 



For the SNLS data set, the available statistical uncertainty was 
used to weight each point so that the equivalent width was cal- 
culated as 

£H^SNLS - (sf , (l - ^] A^,w,] — (3) 

I '='1 7S(W,) 

where w, is the weighting function at each point: 

w,- = l/fr^(^,)' (4) 

The statistical errors were calculated as defined in Equation |2l 
The eiTor induced by any variance in the pseudo-continuum was 
computed with the method described in Sect. 13.11 A description 
of the eiTors caused by sky noise and host galaxy contamination 
is presented in the next sections. 

Note that if either the lower or upper range of an EW mea- 
surement (defined in TableO was outside of the observed wave- 
length range then the measurement had to be removed from fur- 
ther analysis. This constraint affected results on a measurement- 
by-measurement basis and caused the removal of 61 EW mea- 
surements in the SNLS data set. This left a total of 217 measure- 
ments from 73 different high-z SNe in the final set of EW results. 
A number of these measurements were later removed for host 
contamination considerations, however (see Sect. 13.2.21 ). 
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3400 3600 3800 4000 4200 I'lOO 

Rest Wavelength (Angstroms) 
SNI.S SN: OSDlaz 

Fig. 3. An example of an EW measurement (for EW{CaII) in 
this case) on a high-z SNLS object. The median filtered spec- 
trum is plotted on top of the observed spectrum. The local max- 
ima of the feature have been found within this median filter. The 
EW is then computed by numerically integrating the weighted, 
un-filtered flux underneath the pseudo-continuum defined by the 
two maxima. 



3.2.2. High-z EW Systematic Errors II - Host Galaxy 
Contamination 

In nearby objects, the angular separation between the host and 
SN candidate is enough that most observations are not aff'ected 
by host galaxy light. This is not the case with distant SNe. 
In previous high redshift spectroscopic studies, this error has 
been addressed by s ubtracting host galaxy templates during 
the A-^-fitting process (iHoweU & Wangll2002l: iHook et a/.ll2005t 
[Garavini et al. 2007|, during the extraction of the SN spectrum 
dBlondin et aZ.ll200 6). or not at all. In reality, the host contami- 
nation tends to 'wash out' the absorption troughs in SNe spec- 
tra, which causes the measured strength of these features to be 
lower than what would be seen in a spectrum free of host galaxy 
flux. The actual effects of the contamination are dependent on 
the redshift and host galaxy type of each object in addition to 
the amount of contamination. The rolling search method of the 
SNLS provides extensive photometric data that was used to es- 
timate the amount of host galaxy contamination at the epoch of 
spectroscopy for each SN. 

Each object in the SNLS sample has epochs of non-detection 
before and after the explosion. These are used as a baseline from 
which to calculate the precentage increase of the SN light over 
that of the host galaxy on each date when photometry was taken. 
A LIO" (diameter) circular aperture centred on the SN is used 
as it is the basis for the percent contamination estimates that the 
SNLS photometric pipeline provides for each SN candidate. By 
interpolating these measurements to the date of spectroscopy, the 
amount of host contamination in the supernova spectrum can be 
estimated. 

In order to do this, the fact that the spectroscopic aperture 
is not the same size and shape as the photometric aperture must 
be taken into account. The spectroscopic aperture is rectangular, 
with a width of 0.75" and a length set by the extraction win- 
dow used during the spectroscopic data reduction. The amount 
of host light entering the slit was simply scaled from that in the 
LIO" photometric aperture by assuming a spatially flat profile. 
The amount of SN light entering the slit was estimated by mod- 
elling the SN signal as a two-dimensional PSF with a FWHM 
equal to the seeing in the appropriate GMOS acquisition im- 



age. The photometrically determined percent increase can thus 
be converted to a 'percent contamination' of host flux in the ob- 
served spectrum. 

The eff'ects of host contamination on the EW results was de- 
termined by artificially contaminating low-z spectra from several 
different SNe la with template galaxy spectra from E, SO, Sa, Sc, 
and starburst galaxies. In each run of the simulation, the range 
of the observed CFHT i'-band was shifted, based the redshift of 
each SNLS object, to where it would fall in the rest frame. The 
i'-band was chosen because this is the wavelength region where 
SNe la are brightest for the redshift range of this survey; most of 
the SNLS observations are thus made in this band. The flux of a 
chosen host galaxy template was then scaled to match the level 
host contamination (estimated as described above) and added to 
the low-z supernova spectrum. The EW was re-measured to give 
the percent error caused by the specific host type and contamina- 
tion amount for every EW measurement on each SNLS object. 

The results of these simulations were then used to correct 
for this systematic error in the high-z spectra. The median value 
of the percent error estimates among the different low-z spec- 
tra was used as the best estimate for the contamination effects 
on each SNLS supernova, EW measurement, and galaxy type 
in question. The uncertainty of this percent error estimate was 
quantified with the standard deviation of the contamination ef- 
fects from different host galaxy and SNe templates. The effects 
within each general host galaxy type were similar, so this un- 
certainty was only 5% - 7% of the measured EW. Including the 
excess statistical uncertainty from reduced S/N, this error means 
that the uncertainties in the high-z EW results are roughly 2-3 
times those in the corresponding low-z set. The galaxy type for 
the SNe was selected based on the observed host features in each 
affected SNLS spectrum. The results of these corrections on the 
measured EW values are displayed in Fig.|5] 

The results of these simulations helped determine that ob- 
jects with greater than 65% host contamination were too biased 
to be included in the final analysis of this data set. This is be- 
cause at this contamination level, a majority of the high-z EW 
results require adjustments greater than 100% of the measured 
EW and this amount of host flux severely alters the shape of the 
SNe spectra so that the local maxima for each feature cannot be 
identified in a manner consistent with the rest of the sample. This 
removed measurements from 19 objects (in addition to the mea- 
surements removed for observed-wavelength restrictions). The 
final set of high-z results thus includes a total 162 EW measure- 
ments from 54 SNLS SNe la. 



3.3. Ejection Velocities 

The distribution of the ions in the homologously expanding SN 
photosphere and the kinetic energy of the thermonuclear explo- 
sion are expressed in the ejection velocities measured in the SN 
spectrum. Measurements of the ejection velocities of the features 
attributed to Call H&K, SII (near ^5640) and Sill (near ^6150) 
— our EW features 1, 5, and 7 , respectively — have previously 
been studied in both nearby ("Branc h etal.\ 11988 : IWells et al. 
1994; Fisher effl/. 1995; Mazzali eTal.\ Il998t Benetti et al. 



2005; MazzaH et al. 2005) and high redshift dHookef a/."2005t 
Blondin et al. 2006) SNe la populations. The redshifts of the dis- 
tant objects that make up the bulk of the SNLS sample are high 
enough that the SII 'W' (/15640) and Sill /16150 features are gen- 
erally not observed, so only the Call H& K feature was studied. 

The ejection velocity, Vej, was calculated as the blueshift of 
the minimum of the feature away from an estimated rest wave- 
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length of 3945.0 A (a weighted average of the two components 
of the Call H&K feature). The contribution of other ions to the 
absorption in the wavelength range of Call H&K (most notably 
Sill), as well the intrinsic distribution of Call across velocity 
space, can give this feature a 'kinked' shape with two possible 
minima. This shape makes it difficult to determine where the 
mininum of the feature is located. This difficulty is compounded 
by the low S/N of the high-z objects. To make a consistent mea- 
surement between the two sets, the entire Call H&K feature was 
fit with a Gaussian function. The best fit was estimated with ro- 
bust, non-linear least-squares mpfitO routine in IDL. The fit 
was weighted to the most dominant minima of the absorption re- 
gion by making small, initial Gaussian fits to the deepest, bluest, 
part of the feature. The wavelength range from this initial fit was 
then expanded iteratively to span the entire absorption trough. 



and the minimum of the final fit was used to calculate Vej as 



('■^ + iy-i 



H- 1)2 + 1 



(5) 



where A„ is the measured featured minimum, Aq is 3945.0 A 
and c is the speed of light. The statistical error for this fit was 
computed by propagating the uncertanties in the final Gaussian 
parameters into Eq. |5] An additional source of error, caused by 
uncertainty in the selection of points used to fit the Gaussian 
curves, was estimated by taking the difference between the cal- 
culated velocities from the initial and final fits. 

Without having measurement uncertainties available for the 
low-z spectra, the statistical errors were not calculated for these 
objects. For the high-z data, the lower S/N meant these spectra 
had to be smoothed before the Gaussian fits to the Call H&K fea- 
ture could be made. Again a 60 A median filter was used. The 
filtered flux was fit with a Gaussian model that was weighted 
by the error spectra. The statistical error in the final Vq values 
was estimated by propagating the estimated redshift uncertain- 
ties through the velocity calculation. 

No corrections were necessary to account for the systematic 
error from host galaxy contamination as this error was minimal 
for Vej measurements made on the SNe within the 65% contami- 
nation limit discussed in Sect. 13.2.21 

4. Spectroscopic Analysis Results 

4.1. Low-z Results 

The results of the low-z EW measurements, plotted against 
rest-frame epoch (relative to B' maximum light) are shown 
in the first 3 panels in Fig. |4] The mean trend (and Icr 
deviation) exhibited by normal SNe la for each feature is 
shown in the grey contour in each panel. This trend was 
derived by making EW measurements on a template core 
normal Type la SN spec trum made from a combin a tion o f 
core normal spectra by iNugent. Kim. & Perlmutterl (|2002|) . 
The mean trends for the 1991T-like and 1991bg-like objects 
were d erived from similarly constructed template Type l a 
spectra ( Nugent, Kim. & Perlmutted (l2002h : IStern et al.\ (E0041) : 
iLevan et ai. (.2005 ). Nugent private communication) and are 
shown where appropriate. The Icr deviation shown on these 
trends was estimated by expanding the error bars on each trend 
until a fit to the appropriate low-z SNe measurements produced 
a reduced ;(f2 of 1.0. This was done as other methods of estimat- 
ing the Icr range would be more dependent on the varying epoch 



see http : / /cow . physics .wise. edu/ ~craigm/ idl/ 



and object sampling of the low-z data set. A constant scatter as a 
function of epoch was thus assumed, hence the lack of variance 
over epoch in Figs.|4]-|6] 

The core normal objects exhibit fairly homogenous behavior 
in these EW plots. Not surprisingly, this parameter space also 
clearly differentiates between the SNe la sub-types. The lack of 
absorption from IMF (intermediate mass elements) in 1991T- 
like objects (which is due to inferred higher temperatures dur- 
ing thermonuclear runaway in these objects) is clearly seen in 
the EW{CaII H&K) and EWjSill) results. Similarly, the addi- 
tional absorption from (primarily) Till in cooler, 1991bg-like ob- 
jects results in a significant difference with regards to EW{MgII} 
measurements. These results are in agreement with the qual- 
itative and quantitative differences that have already been es- 
tablished for these di fferent S Ne sub-ty pe 

1992a b; Benetti ef a/.l l2005; Bra nch et a/.lB006 

The changes of Call H&K Vq with rest-frame epoch is shown 
in the fourth panel of Fig. |4] The mean trend was estimated by 
fitting a power law to all of the data points from core normal 
objects. The core normal objects display a broad trend about 
this best fit and the SN sub-types are generally interspersed 
around it. Previous estimations of SN la ejection velocities im- 
plied tha t over-luminou s objects have significantly higher veloc- 
ities (Jef frey et a/.ll992|). bu t the statistical analysis of Sill veloc- 
ities by Benetti et al.\ ( 2005b suggests that these objects perhaps 
differ in their velocity gradients, which were not studied here as 
corresponding measurements are not available on the high-z set. 

4.2. High-z Results and Comparison to Low-z SNe 

The EW measurements on the high-z data before and after cor- 
rection for host galaxy contamination are shown in Fig. |5] The 
corrected results are compared to the corresponding low-z results 
in Fig. |6] and all of these high-z results are listed in Table lAl in 
the Appendix. 

Figure |6] shows that the general trends exhibited by the low- 
z sample are indeed followed by the high redshift objects and 
all of the SNLS results are within the range of the low-z SNe 
measurements. A few spectra (marked with open circles in Fig. 
|6]l had EW results and other characteristics that are more similar 
to the SN 1991T-like outliers; these objects are discussed in more 
detail in Sect.|43] 

The fraternity between the high redshift SNe measurements 
and the trends exhibited by the local, core normal objects was 
further investigated with tests. This test compared the high- 
z measurements to the core normal low redshift 'model' illus- 
trated for each feature in Figs. |4] - 16] When making these 
calculations, the full uncertainties for each measurement were 
used. The Icr dispersion in the low-z trends was also considered 
by adding this error (in quadrature) to the uncertainty of each 
high-z measurement. The uncertainty in the estimated epoch for 
every high-z result was also added (again, in quadrature) accord- 
ing to the low-z trend for each measurement. The results, shown 
in Table |4] indicate that there are no statistically significant dif- 
ferences between the spectroscopic features of these SNe la for 
three of the parameters studied here. 

The exception to this result is the comparison of EW{MgII} 
measurements; however, the reality of this difference is ques- 
tionable. It is possible that this difference is a statistical artifact 
caused by the differences in the epoch sampling and total num- 
ber of objects between the high-z and low-z data sets. Figure |4] 
shows that there are an appreciable number of low-z core normal 
points between and -I-IO days past maximum (where the bulk 
of the outliers in the SNLS set are located) which define the Icr 
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Fig. 4. EW Measurements for the 'Call H&K' (top, left), 'Sill' (top, right), and 'Mgll' (bottom, left) features and Call H&K 
Vej (bottom, right) plotted versus rest-frame epoch, for the low-z sample. The shaded, grey region marks the mean trend and Icr 
distribution of the epoch evolution of these measurements in core normal objects. The striped region marks the same trend for 
1991T-like objects in the EW{CaII} and EWfSill} plots and the trend for 1991bg-like objects in the EW{MgII} plot. These trends 
were estimated with EW measurements on template SNe la spectra from the literature. 



Table 4. A quantitative comparison of the spectroscopic features 
in high-z supernovae (as measured with EW and Call H&K ejec- 
tion velocities) and the mean trends exhibited for these quantities 
by local SNe la. These results show that there are no significant 
differences between distant SNLS objects and nearby SNe, with 
the exception of EW{MgII}. 



Feature 


Reduced" 


DOE" 


Prob"^ 


EWICall) 


1.10 


30 


0.67 


EW{SiII) 


0.47 


40 


< 0.002 


EWIMglll 


1.86 


44 


>0.999 


Vej 


1.03 


45 


0.64 



a — X' Idof calculated for the measurement/data set in question and the 
corresponding low-z mean trend 

b — degrees of freedom used when calculating the reduced x'^ 

^datapoints ~ 1 

c — probability (0.00 - 1.00) that the null hypothesis, that the high- 
z data are consistent with the trend displayed by low-z SNe, can be 
rejected. 

range on the empirically derived mean trend. Yet these points are 
from 15 different SNe only; half of the number of SNLS objects 
in the same range. This sampling difference, combined with the 
fact that the low-z sample is not well defined (i.e. it is not a com- 
plete sample) could mean that the true amount of variance in this 
feature is underestimated in the low-z mean trends, which would 



artificially increase the significance of the measured difference 
in the high-z results. If this epoch selection effect is not respon- 
sible for the difference then it is interesting to consider possible 
physical interpretations, which we discuss in Sect.|5] 

4.3. High-z Results: Identification of SNe la Sub-Types 

The comparisons in Fig. |6] show that a handful of the high-z 
SNLS SNe la exhibit spectroscopic features that are more con- 
sistent with the 1991T or 1991bg-like SNe from the low-z sur- 
veys. 

Overluminous 1991T-like SNe la have less absorption from 
IME at certain epochs. Three outliers in the EW{Call) and 
EW{Sill} results have evidence of being similar to the known 
class of 1991T-like objects. Their spectra are shown in Fig.|2] 
03D4cj displayed a smaller amount of Sill than a typical core 
normal object. The error bars on this measurement are relatively 
large (EW{Sill}= 5.9 + 6.1 A), but the photometric data for 
this object support s the conclusion that it is a 1991T-like SN; 
lAstier et al. I ([2006') stated that the lightcurve from this SN was 
incompatible with the model used to fit all of the other SNLS 
SNe. The;t^^ template fits to 03D4cj also indicate that it is simi- 
lar to SN 1991T. 

The 199 IT- like classification also appears valid for SNLS 
supernova 03Dlcm. This SN has a very low EWjCall) mea- 
surement and an EWjSill) result that (within large error bars. 



10 



T. J. Bronder et al.: SNLS Spectroscopy: Testing for Evolution in Type la Supemovae 




Fig. 6. The results from the SNLS high-z SNe (open squares) compared to the EW results from the low-z sample (filled symbols). 
The contours in each plot are the same as Fig.|4] The objects that were identified as 1991T-like SNe are circled in red. 



EW{SiII) = 14.8 + 10.6 A) is also consistent with the low val- 
ues measured in 1991T-like objects. Photometrically, 03Dlcm 
is also more like a 1991T-like SN as it had one of the high- 
est stretch values f rom the entire set of Gemini SNLS data at 
s = 1.173 + 0.061 dAstier et al. Il2006h . The best matching tem- 
plate spectrum also corresponds to an overluminous low-z SNe. 
The spectral and photometric dates do differ (+1 day and -4.42 
days, respectively) in this case, but the low S/N on this object and 
inherent uncertainty on spectral dating with this template match- 
ing method (which is on the order of a few days) means that this 
is not necessarily enough evidence to outweigh the 1991T-like 
EW and s values of OBDlcm. 

The SNLS SN 04D3mk also had a very low, well constrained 
EW{SiII) value. At the epoch this object was observed, the 
EW{Mgll) result cannot differentiate between normal and over- 
luminous SNe and the EW{Call} measurement was unusable as 
it fell on the gap between GMOS detectors (e.g.. Sect. 13.2b . The 
final photometry of this object did not reveal any photometric 
peculiarities, but the x' template matches to 04D3mk singled 
out overluminous SNe la spectra only as the most similar low- 
z objects. This supports the quantitative similarities of 04D3mk 
to this particular SN la sub-type. With this evidence, SNLS ob- 
jects 03Dlcm and 04D3mk are the most distant supernova (at 
z = 0.87 and z - 0.812, respectively) to be shown to be spectro- 
scopically similar to the 1991T-like SNe la sub-types. To reiter- 
ate, the primary source of this identification for these objects is 
the EW measurements which are quantitatively more similar to 
the overluminous sub-class of SNe la than other sub-types. 

The SNLS SNe 03Dlco and 05D3mq are 1991bg-Hke can- 
didates as their large EW{MgIl) values are more consistent with 



the excess absorption measured in the subluminous low-z SNe. 
However, these two high-z SNe had measurable amounts of Sill, 
a feature that has yet to be observed in any 1991bg-like SNe be- 
cause it is obscured by additional absorption from Till and Fell. 
This evidence, combine d with the normal photometric param- 
eters (A stier et al. II2006I) derived for these objects and the fact 
that template matches to 1991bg-like objects resulted in poor;^^^ 
results, is enough to justify that these SNe are more like core nor- 
mal SNe than any of the underluminous SNe observed to date. 

4.4. Correlations With Host Galaxy Type 
4.4.1. Low-z spectra 

The general trend of brighter SNe and a higher dispersion of 
SNe la peak magnitudes in late-type galaxies is w ell documented 
(iHamu v efaZ.II 19961 l2000t ISulhvan et a /."2003') . Yet only a few 
quantit ative spectroscopic corollaries (Branch & van den Berghl 
il993l:IBranch. Romanishin. & Ba ron 1996) have been made that 
agree or disagree with this trend. Kolmogorov-Smirnov (KS) 
tests were thus used to explore any differences between the spec- 
troscopic parameters measured in SNe from different host types 
in our low-z sample. This test was carried out on the residuals 
of these measurements, where the residual for each EW and Vej 
value is defined as the difference between the measurement and 
the corresponding low-z mean trend. 

The host type for each low-z SN was grouped into the broad 
categories of 'E/SO' and 'Spiral' (Sa, Sb, Sc and irregular) 
galaxies. The published morphology (from the Nearby Galaxy 
Catalogue and other appropriate databases) of each SN host 
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Fig. 5. Measurements of the EW of the 'Call H&K', 'SilF, and 
'Mgir features (from top to bottom), versus rest-frame epoch, 
in the SNLS high-z sample. The shaded and striped regions are 
the same as in Fig. |4] The EW measurements that did not need 
any correction for the systematic effects of host galaxy contam- 
ination are shown with the dark, filled squares. All of the other 
points illustrate the measurements that were corrected for host 
contamination. The adjusted EW values are connected to the cor- 
responding initial measurements with dotted lines. 



Fig. 7. Three high-z SNe la - 03D4cj, 03Dlcm and 04D3mk - 
that demonstrate a quantitative similarity to the 1991T-like SNe 
la sub-type. The green spectrum is the observed spectrum with- 
out any binning. The black spectrum is the best-fitting SN from 
the low-z library and the red spectrum a smoothed version of the 
observed spectrum. The primary evidence for the similarity of 
these objects to 1991T-like SNe is the EWfCall) and EWjSill) 
measurements, although in each case the fits to SN 199 IT 
or SN 1999aa were quantitatively better than those of core nor- 
mal SNe la. Photometric information also supports this sub-type 
identification for these three SNe. 



galaxy (based on the hosts listed in the SNe reference papers, 
see Table[T]) was used to determine which group was appropriate 
for each SN host. 

A comparison of the residuals for the four measurements dis- 
cussed in this paper, grouped according to these two broad host 
galaxy groups, are displayed in Fig. [8] At low-z the only feature 
that shows a significant trend with host type is Sill. 



4.4.2. High-z spectra 

Another test for possible evolution is to determine whether the 
correlations between low-z SNe la characteristics and host mor- 
phology are observed within the high-z set. The host types were 
estimated for the SNLS objects that had clear host galaxy fea- 
tures in their spectra and placed into one of three broad cate- 
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Fig. 8. Histograms comparing the residuals of EW and Call Vej measurements in SNe from E/SO (open histograms) and spiral 
(filled histograms) host galaxies. The residuals from low-z SNe are displayed in the top row and the SNLS high-z results are on 
the bottom. The vertical dashed line in each panel is the mean value for each distribution. The histograms have been normalized to 
have the same area. The distribution of low-z SNe results show a significant difference for the objects from early and late hosts (as 
determined with KS tests). The high-z EW{SiIl) results show a similar trend while EW{Mgll) distributions display a difference that 
is not seen at low-z. 



Table 5. Comparison of SNe la spectral features in different 
host types at low and high-z. The values given are the mean and 
error in the mean of the residuals from the low-z trend. 



Feature 


SNe Measurement (low-z trend subtracted) 




E/SO 


Spiral 


Low-z 


EWICall) - A 


-1.7 ±4.7 


-2.0 ± 3.2 


EWISilll - A 


5.1 ±0.7 


-4.5 ± 0.9 


EW{MgII) - A 


-0.6 ±3.8 


-0.7 ± 2.9 


Vej - 1000 km/s 


-0.18 ±0.23 


0.17 ±0.27 


High-z 


EWICall) - A 


-13.0 ± 16.0 


-0.8 ±7.2 


EWISill) - A 


5.0 ±2.4 


-1.4± 1.1 


EWIMgll) - A 


29.1 ±9.5 


0.1 ±7.7 


Vej - 1000 km/s 


-0.13 ±0.90 


0.55 ± 0.50 



gories; 'E/SO' for early (elliptical and SO hosts), 'Early Spirals' 
for objects with spectroscopic features similar to Sa and Sb 
galaxies, and 'Late Spirals' for Sc and irregular hosts. The sub- 
set of SNLS SNe la with indentifyable host galaxies includes 45 
of the 54 SNLS objects that had EW or Vej results. 

The distribution of these SNLS measurements, displayed in 
the histograms in the bottom of Fig. [8] do show clear similari- 
ties to the low-z distributions (top, Fig.[8]l. The distant late type 
galaxies host SNe that display a smaller amount of Sill than ob- 
jects from early hosts and the ejection velocities measured in 
spiral galaxies peak at a higher value. A comparison of the mean 
values for each measurement (again after subtracting the low-z 
trend to remove dependence on SN phase) divided by host type 
supports this conclusion, as can be seen from Table |5] 



Similar to the low-z results, the high-z data show a trend with 
regards to host type and EWfSill), though at a slightly lower sig- 
nificance because of the larger high-z uncertainties. We also see a 
possible difference in Mgll that is not seen at low-z. Thus, the KS 
tests of the distribution of high-z EW(Sill) and Vej results do not 
show as significant of a statistical difference between SNe from 
different host types. This is probably due to the smaller num- 
ber of measurements in the high-z data in addition to the larger 
errors. For example, the measured maximum deviation between 
the EW{SiIl} distribution in the high-z galaxies (KS statistic D) 
is D - 0.50. This value is comparable to the low-z D = 0.64 and 
would indicate a significant difference if the high-z sample was 
larger by only 10 measurements. 

The trend for lower EWs and higher ejection velocities in 
late type galaxies seen at low and now also high-z is a corrolary 
to the broad-ban d photometric s t udy o f SNLS SNe host galax- 
ies completed by ISullivan et al.\ (l2006l) . which showed that the 
lightcurve shapes of high-z SNe demonstrate a dependence on 
host galaxy star formation rates in a manner consistent with low- 
z objects. 

4.4.3. Dependence on Redshift 

The high-z EW and Vej measurements (after subtracting the low- 
z trend to remove the dependence on SN phase) are plotted 
against redshift in the left-hand column in Fig. |9] No depen- 
dence on redshift is seen. Calculations of Spearman's p rank 
correlation and Kolmogorov-Smirnov tests attest that there are 
no trends in the data with redshift that are significantly different 
from random scatter. Thus we see no evidence in our sample for 
systematic evolution of SNe la properties in the large redshift 
range covered by the SNLS. 
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Fig. 9. Left: The high-z residuals (i.e. the difference from low-z SNe la core normal mean trend) from the EW and Vgj measurements 
are plotted against redshift. Right: the same data divided by host galaxy type. 
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Fig. 10. A comparison of the EW measurement of the 'Sill 
4000' feature to the absolute B' peak magnitudes of low-z SNe 
la (black circles) and distant SNe (lighter circles). The signif- 
icant correlation in the low-z data reveals a new spectroscopic 
sequence in SNe la that can be used to estimate SNe la peak 
magnitudes. The high-z points have larger error bars with re- 
gards to both EW and magnitude, but are still consistent with 
the low-z trend (illustrated with the solid line). 

On the right side of Fig. |9l the measurements are divided 
by host galaxy type. The evolutionary histories of host galaxies 
within a given morphological category have more in common 
than galaxies of different classfications. The comparison of SNe 
la spectroscopic features at different redshifts but within similar 
host types should thus be more sensitive to any evolutionary ef- 
fects than tests using the bulk properties of the high-z sample as 
a whole. The panels to the right in Fig. |9] illustrate that there are 
no redshift-dependent trends within the SNe results from each 
host type. 

4.5. EW{Sill] - A New SNe la Spectroscopic Sequence 

During the course of this study, an interesting correlation was 
found between EW{Sill) measurements and Mb for the low-z 
SNe. Figure [TO] shows the low-z SNe EW{Sill) measurements, 
taken as close as possible to maximum light (within +7.0 days), 
compared to the absolute peak (B') magnitudes (uncorrected 
for any lightcurve shape-luminosity relationship) that were esti- 
mated for 13 of the low-z SNe (see Sect. l2.2.2] i. A Spearman rank 
correlation test found a correlation with 3cr significance between 
the measured EW{Sill) and calculated peak magnitude, Mb. The 
best fit linear correlation of these quantities is described by the 
following: 

Mb = (-20.07 + 0.19) + (0.039 + 0.01 1) x (EWfSill}). (6) 



spectrum out to redshifts of about 1 .2 [unlike measurements in 
the redder end of the spectrum such as the A5912 feature studied 
in|Hachinger, Mazzali, & Benetti (2006)]. 

One major physical parameter that underlies this se- 
quence is tempera ture (Nugent et a/.lll995at [Branch et a/.ll2006l: 
iBaron et'al. 1 120061) . For an assumed Chandrasekhar mass (Mch) 
WD explosion scenario, the distribution of peak magnitudes is 
based primarily on the different amounts of ^''Ni synthesised. 
Higher production of this isotope, which is made at the expense 
of IMF such as Sill, results in a brighter SNe and causes the 
ejecta to be hotter. This tends to overionise the already smaller 
amount of IMF observed at maximum light thus results in the 
appearance of a smaller amount of Sill in the spectra of hotter 
and brighter SNe at these epochs. 

Evidence that this spectroscopic sequence is similar at high- 
z is illustrated in the SNLS points in Fig. (TO] In this Fig. the 
absolute magnitudes of the SNLS high-z SNe, uncorrected for 
lightcurve stretch or colour adjustments, were computed as- 
suming the best fit cosmology and published values from 
lAstierefg/.l (120061) . 

Within the relatively large EW error bars, the SNLS super- 
novae are still consistent with the low-z trends in this parameter 
space. The next logical step is to determine whether the appli- 
cation of this relation to high-z SNe, just like the use of s or 
Am,, parameters, decreases the scatter of luminosity distances in 
a Hubble diagram. This test is displayed in Fig.fTTI 

For this test, the set of SNLS SNe that were published in the 
first year SNLS cosmology results and had measured EW{SiII) 
values were used. A total of 18 SNe la from lAstier et al. 1 (120061) 
met these criteria. The high-z EWjSill) measurements were used 
with the relation estimated from the low-z sample (Eq.|6]l to esti- 
mate the absolute peak B-band magnitudes (now termed Msiu). 

The distance modulus (//) for these distant objects was then 
calculated as 

ju = niB - Msiii, (7) 

where niB is the o bserved peak B' ma gnitude previously pub- 
lished in Table 9 o f lAstierefg/.l ( 12006 1). The uncertainties in the 
linear MB-EW{SiII} fit and niB were propagated along with the 
EW errors during this calculation. 

Figure[TT]indicates that the EW{Sill) sequence does provide 
a plausible method to use for estimating peak SNe la distances 
for cosmology. The dispersion in the EW{Sill}-estimat ed dis- 
t ances about the best-fit cosmological model from Asti er et al. I 
( |2006|) is 0.26 mag., which is near the computed residuals of 0.23 
mag. from the lAstier et al. I (12006) stretch and colour-corrected 
SNLS SNe. Both of these values are much improved over the 
0.37 magnitude dispersion from the un-corrected SNe magni- 
tudes in this sub-set of 18 high-z SNLS SNe. The application 
of this correlat ion to SNe la cosmology is further discussed in 
iBrondeil ( 120061) . 



The small residuals from this fit place this spec- 

tral indicat or o n par with others in the literature 

dwells et al.\ Il994t iBongard etal.\ l2006t iNugent ef a/.l Il995at 
iHachinger. Mazzafi. & Benettil l2006h . This spectroscopic se- 
quence is a lso relatively free of the epoch const raints on similar 
sequences dNugent et al.\ Il995al; iBenetti et al] 120051: IKMI) be- 
cause EW{SiII) evolves little near maximum light and only one 
measurement is needed to constrain Mb with this correlation. 
Any single EW measurement near this date (the authors found 
that +7.0 days was best as the EW{Sill) values evolved little 
over this period) yielded the similar results. The wavelength 
of this feature also makes it ideal as it falls within the optical 



5. Discussion 

As discussed in Sec. |4l the spectroscopic measurements and 
comparisons completed on this new set of distant SNe have 
shown that these SNLS-observed objects display an overall sim- 
ilarity with the SNe la observed at low redshifts. We find no 
significant difference between the SNe la samples for three of 
the four features studied and a possible difference in the fourth 
feature (EW{Mgll)). We also find similar trends with host type 
at low and high redshift. Below we discuss these results in the 
context of the population of progenitors in addition to the rates 
of various Type la sub-types in the SNLS high-z sample. 
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Fig. 11. Hubble diagrams for a subset of SNe la from the SNLS first year dataset. Solid points show the distance moduli uncorrected 
for stretch or colour effects. Left: The lighter (or coloured) points show distance moduli that were calculated with peak magnitudes 
estimated by the EW{SiII) sequence. Right: The lighter (or coloured) points show distance moduli that were corrected with s (stretch) 
and c (colour) corrections ( Astier et al. 2006). The solid line indicates the best fit cosmology from the first year SNLS results - 
a spatially flat Universe with = 0.263 and Qa = 0.737. The dashed line is the model from a flat Universe with Dm = 100. 
The residuals about this cosmology are shown in the bottom two panels; the uncorrected high-z SNe points have a dispersion about 
this cosmology of 0.37 mag nitudes, whil e the E WjSilll points have a dispersion of 0.26 magnitudes; this is comparable to the 0.23 
magnitude dispersion in the lAstier et al. I ( l2006i) results. 



5.1. Low-z and High-z Discrepancy - EW{Mgll} 



In Sect. 14.21 a possible difference in EW{MgII} between low 
and high-z samples was presented. For the reason given in that 
section (namely the differences in the time-sampling of obser- 
vations in the two samples), the significance of this difference 
may n ot be hig h. This was also sugested in a similar study by 
iGaravini et a/.l (l2007l) . who also noted the need for an analysis 
of a larger sample of normal SNe la to more accurately deter- 
mine the true mean trend for this measurement. 

However, it is interesting to consider possible physical 
explanations, for example it is possible that it is a sign 
of metallicity differences. Explosion and nucleosynthesis cal- 
culations from iTravaglio. Hillebrandt. & Reineckd (l2005h and 



iRoepke et a/.l (l2006 l). which were computed by modelling ther- 
monuclear burning by deflagration (as opposed to delayed deto- 
nation or detonation only) of a C-hO WD progenitor, showed that 
Type la SNe from higher metallicity (Z) environments produce 
a significantly lower (by a factor of 10) amount of ^"^Mg. In their 
results, the abundances of Si and Ca were relatively unaffected 
by metallicity changes. Thus, it may be that the group of high-z 
SNe that are below the mean trend for this feature are all SNe 
from higher metallicity progenitors. The normal EWjSill) and 



EW{CaII} measurements in these SNe support this explanation 
within the context of these models. 

Conversely, model S Ne la spec t ra fro m progenitors of 
varying Z calculated by iLentz et al] (l2000l) showed the op- 
posite eff'ect in the EW{MgII} measurements. The iLentz etaU 
(l2000h spectra were calculated with a similar deflagratio n 
model SN [W7 from iNomoto. Thielemann. & Yokoil (1 1984 1. 
with the metallicity changes approximated by scaling the num- 
ber abundances of elements heavier than O in the unburned 
C-hO layer. The fi nal IME abundances were not published in 
iLen tz et aP (lOOO h, so they could not be com pared directly to 
the Travaglio, Hill ebrandt. & Reineckd (l2005h results. Instead, 
the EW measurements defined for this study were made directly 
on the final Lentz et al. (2000) model spectra. These measure- 
ments showed that the SNe with higher Z had larger EW{MgII} 
measurements at pre-maximum epochs but then evolved such 
that the diff'erent Z models were nearly indistinguishable after 
maximum light. 

Another factor that may be contributing to this difference is 
the tr end for the SNLS to find more slightly bluer, brighter ob- 
jects (As tier et ani2006 ). F04 noted that the evolution of diff'er- 
ent SNe along the mean trend for EW{MgIIj correlated with lu- 
minosity, with brighter objects making the 'jump' to higher EW 
at later times. The small group of outliers below the EW{MgII} 
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mean trend between +5 and +10 days may just be reflecting a 
slightly brighter group of SNe in the high-z sample. The average 
s value within these SNe is just above 1 .0 [th is corresponds to 
brighter SNe, see Fig. 9 in lAstier et al. I (l2006l) 1. which supports 
this hypothesis. This feature comparison may be more sensitive 
to slight differences in the strecth distribution of the samples than 
the other EW results. Further investigation into this parameter, 
both observationally and theoretically, is certainly warranted. 

5.2. Rate of Peculiar SNe 

The EW results of this large high-z set note 2 to 3 objects (out of 
54 distant supernovae with EW results, ~ 3 - 6%) that are spec- 
troscopically similar to the class of 1991T-like SNe sub-types 
that have been observed locally. Previous studies of distant SNe 
la have noted a similarly small percentage of these sub-types in 
the high-z samples; only spectra from SN e 2002iv (z - 0.23 1 ) , 
2003jn (z = 0.33), 20 02fd (z = 0.27) jGaravini ef a/."2007|) 
and 2 0031h (z = 0.56) dMatheson et a/.l l2005; Krisciunas et alj 
12005 ) have been shown to be quantitatively similar to 199 IT or 
1999aa. This is in contrast to the apparent contribution of these 
sub-t ypes to the low-z SNe sample, whic h is roughly 10% to 
20% dLi et a/.ll200'Ti IStrolger et a/.l l2002l). This has be en inter- 
preted as a possible sign of evolution jBranchlbOOU iLi et al.\ 
1200 la'b') or of a significant difference between the nearby and 
distant SNe la populations. The identification of the 3 possible 
1991T-like objects in this high-z sample nearly doubles the num- 
ber of these sub-types observed at higher redshifts, but the per- 
centage of identified 199 IT- like SNe is still much lower than at 
low-z. It is difficult to determine the significance of this apparent 
discrepancy due to the unknown influence of selection effects. 
For example, the low-z peculiar objects are often singled out for 
spectroscopic follow up and thus may be over-represented in the 
nearby sample. Another complication is the uncertainty with re- 
gards to whether the overluminous 1991T-like and 1999aa-like 
SNe are distinct sub-classes of Type la SNe or if they represent a 
possible continuum of thermonuclear SNe dBranch et al.WlOOdti . 
Note that the EW and Vej measurements cannot readily address 
this latter issue; the overluminous high-z objects are referred to 
as '1991T-like' in the interest of brevity only and does not rep- 
resent an attempt to clarify the debate on the 1991T-1999aa re- 
lationship. 

It may be that the difference between the fraction of these 
sub-types at different redshifts is primarily an effect of the epoch 
coverage differences between these data sets. SNe 1991T and 
1999aa displayed unusual spectroscopic features primarily be- 
fore maximum light and then evolved to look little different than 
core normal SNe. Within the low-z set considered here, 77% of 
the SNe were observed at least once prior to maximum light, 
compared to only 33% (18 of 54 SNe) of the high-z sample. With 
the simple assumption that this factor of ~ 2 difference has re- 
duced the number of overluminous SNe possibly identified in the 
high-z sample by a similar amount, the disparity in the percent- 
age of overluminous sub-type identifications between the two 
data sets is explained. 

Even if this explanation is true, there is still a lack of of 
the underluminous 1991bg-like SNe at high redshifts. This is 
perhaps expected for this data set as many factors work against 
SNLS spectroscopic identification of these objects. Mostly, it is 
because they are intrinsically much fainter and thus harder to 
observe photometrically and spectroscopically. This implies that 
even if the SNLS is observing these objects spectroscopically, 
they would be much harder to conclusively type. Thus the lack 
of 1991bg-like SNe la in this set is not unexpected. 



The end of Sec. 12. 1 . 1] summarizes that, due to the SNLS sur- 
vey setup, high-z SNe with vastly different properties than the 
nearby objects would most likely not be selected for spectro- 
scopic observation. Thus, the conclusions of the tests made here 
cannot rule out that such a population of SNe la exists at high-z, 
although none of the results show any indication of new sub- 
classes of distant SNe. However, we note that occasionally very 
rare and unsual obj ects are found [see SNLS object 03D3bb in 
iHowell et aU (l2006h . for example]. 

5.3. Implications for Progenitors 

Comparisons of SNe la across different host galaxy types can 
also reveal keys to the progenit or scenarios for these objects . 
One model of SNe la formation ( Scannapi eco & Bild sten 2005h 
has used the implic ations of o bserved SNe la enviro nmental de- 
pendencie s (Wang et al. 1997; Umeda et al. 1999; H amuv et al.\ 
l2000tlHowell.200UlMannucci et al.,2005!) to create a model that 
describes SNe la as two distinct populations; a 'prompt' com- 
ponent dependent on the galaxy star formation rate and a 'de- 
layed' component whose rate is proportional to galaxy mass. 
This model has been supported with SNLS data (ISulhvan et al.\ 
12006 !) and may be useful in predicting any evolutionary shifts 
in Sne la properties (Howell et al. in prep). The measurements 
made in this study can also be put into the context of this 'two- 
component' model. Tableland Fig|8]show a trend of Sill with 
host type (at both low and high-z) in the direction expected 
from this two-component model. At all redshifts, SNe in spiral 
hosts show weaker Sill. Assuming that these early type galax- 
ies are star-forming, and hence host the 'prompt' SNe la (which 
also tend to be the brighter, higher stretch SNe la), and that 
the EW{Sill)-MB correlation is valid, then the trend of lower 
EW(Sill} measurements in SNe from late type galaxies agrees 
with the main aspects of the two-component model. 

Another implication of this model is that as observations 
are extended to higher redshifts, a larger sampling of SNe la 
from the 'prompt' star formation-based population would be ex- 
pected. As these 'prompt' SNe are brighter than their 'delayed' 
counterparts, this should correspond to additional observations 
of SNe with high s values and correspondingly low EW{Sill) 
measurements. This effect can be tested in the SNLS sample 
considered here by combining the high-z results from all galaxy 
types, tests and measures of correlation show no statistical 
difference in EWjSill) between low and high-z samples (see also 
Fig.|9]l. However, this could be due to a lack of sensitivity to de- 
tect this subtle drift in the population. A possible shift of 10% 
in the brightness (or stretch) of the high-z population (Howell et 
al. in prep) would correspond to a change of the order of a few 
(x! 4) A in EWjSill); Table |5] shows that the error in the means 
of the distribution of high-z EW(Sill) is on the order of this shift 
and is much higher than the corresponding error in the low-z re- 
sults. An SNLS sample with 2-3 times more SNe (which will be 
available before the survey is complete) should be able to detect 
this shift toward the 'prompt' SNe. 



6. Conclusions 

New results from a quantitative study of a large set of spectra 
from high-z SNe la used for cosmological fitting have been com- 
pared to corresponding measurements on a sample of SNe la ob- 
served locally. The measurements from all of the studied high-z 
spectra are within the range of the results at low-z. A statistical 
comparison of these data sets reveals that there is no significant 
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difference (less than 2cr) in the strength of 2 SNe la-specific ab- 
sorption features and the ejection velocity of Call H&K. The 4* 
feature - which quantifies absorption from Mgll and other ions 
near 4300 A - shows a possible difference that should be inves- 
tigated further 

These methods also reveal a new spectroscopic sequence in 
SNe la based on the strength of the Sill feature near 4000 A at 
epochs close to maximum luminosity. An initial test of this se- 
quence was made on a sub-set of 18 SNLS SNe that had the ap- 
propriate EW{Sill) measurem ent and were published with previ- 
ous SNLS cosmology results (lAstier et al. 2006). This demon- 
strated that the EW{Sill)-MBp^,t relation reduces the scatter in 
high-z luminosity distance estimates in a manner consistent with 
the lightcurve shape-luminosity corrections that are currently 
used to standardise SNe la peak magnitudes for cosmology. 

At present this type of spectroscopic analysis is one of the 
most promising avenues of investigating the dispersion, homo- 
geneity and possible evolution of distant Type la supernovae. 
This paper has used empirical techniques on the largest set 
of distant SNe la data studied to date. Although the selection 
method used to generate the SNLS high-z sample implies that 
the existence of a new population of Type la SNe at high-z with 
unusual spectral properties cannot be ruled out, a large set of 
high-z SNe have been shown to be spectroscopically similar to 
nearby objects. At high-z the same trends of SN properties with 
host type are seen as at low-z, and furthermore, within each host 
galaxy class no evidence for evolution in SN properties with red- 
shift is seen. These results represent an important step in ad- 
dressing one of the most pressing issues for Type la cosmology. 
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Table 6. Observed properties and instrument settings for SNLS SNe candidates observed from November 2004 to May 2006. The magnitudes and 

"percentage increase" are estimated from the lightcurves and correspond to the date of spectroscopy. 



SN Name 


RA (2000) 


Dec (2000) 


UT Date 




Exposure (s) 


Mode 


'^cent 


seeing (") 


Mag (i') 


%mc 


04Dlpu 


02:27:28.437 


-04:44:41.71 


2004 


-12 


-12 


7200 


N&S 


720 


0.49 


24.08 


39 


04D2mh 


09:59:45.872 


+02:08:27.94 


2005 


-01 


-09 


2400 


C 


720 


1.00 


22.90 


730 


04D2mj 


10:00:36.535 


+02:34:37.44 


2005 


-01 


-09 


3600 


N&S 


720 


0.62 


23.19 


54 


OSDlaz 


02:25:12.492 


-04:36:08.17 


2005 


-09 


-04 


5400 


N&S 


720 


0.63 


23.78 


1246 


OSDlby 


02:24:35.448 


-04:12:04.16 


2005 


-09 


-03 


5400 


N&S 


720 


0.62 


23.56 


9 


OSDlbz 


02:25:00.685 


-04:12:00.74 


2005 


-09 


-02 


7200 


N&S 


720 


0.88 


23.96 


4310 


05D1CC 


02:26:31.173 


-04:09:53.40 


2005 


-09 


-05 


3600 


N&S 


720 


0.46 


23.18 


114 








2005 


-09 


-06 


1800 






0.83 


23.18 


114 


OSDlee 


02:24:46.766 


-04:00:01.86 


2005 


-10 


-31 


3600 


N&S 


720 


0.59 


23.24 


25 


OSDlej 


02:26:06.317 


-04:43:45.79 


2005 


-10 


-27 


3600 


N&S 


720 


0.47 


22.83 


15 


OSDlem 


02:24:05.503 


-04:56:23.35 


2005 


-11 


-05 


3600 


N&S 


720 


0.62 


23.71 


476 








2005 


-11 


-06 


1800 






0.57 


23.71 


9524 


OSDleo 


02:25:23.116 


-04:38:33.36 


2005 


-10 


-27 


5400 


N&S 


720 


0.39 


23.75 


21 


OSDler 


02:25:41.466 


-04:00:16.32 


2005 


-10 


-28 


7200 


N&S 


720 


0.55 


23.99 


94 


OSDlhn 


02:24:36.254 


-04:10:54.94 


2005 


-12 


-04 


2700 


C 


680 


0.63 


20.86 


249 


OSDlju 


02:26:57.303 


-04:08:02.78 


2005 


-12 


-27 


7200 


N&S 


720 


0.50 


24.04 


91 


OSDlka 


02:24:58.945 


-04:34:28.56 


2005 


-12 


-28 


5400 


N&S 


720 


0.86 


23.21 


1059 


OSDlkl 


02:24:33.544 


-04:19:08.33 


2005 


-12 


-31 


5400 


N&S 


720 


0.59 


23.61 


247 


05D2ab 


10:01:50.833 


+02:06:23.02 


2005 


-01 


-09 


1800 


C 


680 


0.67 


21.83 


263 


05D2ah 


10:01:28.704 


+01:51:46.18 


2005 


-01 


-09 


1800 


c 


680 


0.96 


20.83 


1017 


05D2ck 


10:00:45.203 


+02:34:22.13 


2005 


-02 


-16 


5400 


N&S 


720 


0.65 


23.88 


7592 


05D2ja 


10:00:03.809 


+02:17:36.12 


2005 


-05 


-07 


5400 


c 


680 


0.94 


22.15 


115 


05D2nt 


10:00:58.234 


+02:22:21.62 


2005 


-12 


-26 


5400 


N&S 


720 


0.69 


23.46 


340 


05D2ob 


09:59:00.705 


+01:50:56.62 


2005 


-12 


-28 


5400 


N&S 


720 


0.73 


24.22 


242 


05D3ax 


14:19:17.595 


+52:41:15.07 


2005 


-02 


-11 


5400 


N&S 


720 


0.45 


23.40 


342 


05D3bt 


14:18:35.980 


+52:54:55.95 


2005 


-03 


-06 


3600 


C 


720 


0.84 


23.52 


71 


05D3cf 


14:16:53.369 


+52:20:42.47 


2005 


-03 


-08 


3600 


c 


720 


0.69 


22.96 


94 


05D3ci 


14:21:48.085 


+52:26:43.33 


2005 


-03 


-09 


5400 


N&S 


720 


0.76 


23.03 


70 


05D3cq 


14:18:46.173 


+53:07:55.55 


2005 


-04 


-11 


5400 


N&S 


720 


0.64 


24.10 


156 


05D3CX 


14:21:06.560 


+52:45:01.70 


2005 


-04 


-09 


5400 


N&S 


720 


0.50 


23.52 


65 


05D3jq 


14:21:45.462 


+53:01:47.53 


2005 


-06 


-04 


5400 


N&S 


720 


0.60 


23.36 


30 


05D3ki 


14:21:16.341 


+52:35:42.47 


2005 


-05 


-14 


5400 


N&S 


720 


0.67 


23.68 


49 


05D3km 


14:22:38.298 


+53:04:01.14 


2005 


-06 


-04 


5400 


N&S 


720 


1.0 


23.89 


536 


05D3kp 


14:20:02.952 


+52:16:15.28 


2005 


-06 


-03 


5400 


N&S 


720 


0.81 


23.36 


7848 


05D3kt 


14:19:53.730 


+52:44:34.64 


2005 


-06 


-03 


5400 


N&S 


720 


0.82 


23.53 


141 


05D31b 


14:17:31.775 


+53:10:04.38 


2005 


-06 


-05 


5400 


N&S 


720 


0.51 


23.38 


192 


05D31C 


14:22:22.902 


+52:28:44.11 


2005 


-06 


-06 


3600 


C 


720 


0.63 


22.67 


17 


05D31W 


14:17:44.021 


+53:06:24.98 


2005 


-07 


-01 


7200 


N&S 


720 


0.45 


24.41 


170 


05D3mn 


14:18:45.206 


+52:19:23.56 


2005 


-07 


-12 


7200 


N&S 


720 


0.77 


23.31 


77 


05D3mh 


14:18:59.805 


+52:40:03.50 


2005 


-07 


-10 


7200 


N&S 


720 


0.54 


23.50 


30 


05D3mq 


14:19:00.398 


+52:23:06.81 


2005 


-08 


-03 


1800 


C 


680 


0.71 


21.48 


4748 


05D3mx 


14:22:09.078 


+52:13:09.35 


2005 


-07 


-30 


2700 


c 


680 


0.83 


23.10 


>10000 


05D3ne 


14:21:02.946 


+52:29:43.92 


2005 


-07 


-14 


2700 


c 


680 


0.85 


20.48 


1074 


05D4av 


22:14:10.515 


-17:54:42.67 


2005 


-07 


-10 


3600 


N&S 


720 


0.53 


23.2 


77 


05D4bm 


22:17:04.621 


-17:40:39.45 


2005 


-07 


-11 


3600 


c 


680 


0.55 


22.04 


144 


05D4ca 


22:14:11.350 


-17:48:15.47 


2005 


-08 


-02 


2700 


c 


720 


0.63 


22.24 


156 


05D4cn 


22:13:31.454 


-17:17:19.92 


2005 


-08 


-02 


5400 


N&S 


720 


0.70 


23.30 


82 


05D4dt 


22:14:25.851 


-17:40:16.03 


2005 


-09 


-02 


3600 


N&S 


720 


0.55 


22.62 


20 


05D4dx 


22:13:39.387 


-18:03:20.99 


2005 


-09 


-06 


5400 


N&S 


720 


0.47 


23.14 


39 


05D4dy 


22:15:30.102 


-18:12:55.54 


2005 


-09 


-06 


5400 


N&S 


720 


0.53 


23.87 


684 


05D4fn 


22:17:01.076 


-17:49:25.16 


2005 


-10 


-08 


7200 


N&S 


720 


0.65 


23.86 


4474 


05D4fo 


22:15:20.925 


-17:16:05.25 


2005 


-10 


-07 


3600 


N&S 


720 


0.77 


22.47 


1845 


05D4gw 


22:14:47.365 


-17:31:54.99 


2005 


-11 


-05 


5400 


N&S 


720 


0.46 


23.69 


283 


05D4hn 


22:17:13.545 


-17:54:45.40 


2005 


-12 


-01 


5400 


N&S 


720 


0.50 


23.50 


474 


06D3bz 


14:17:10.044 


+53:01:29.31 


2006 


-02 


-08 


3600 


N&S 


720 


0.71 


23.20 


115 


06D3cb 


14:20:43.536 


+52:11:28.11 


2006 


-02 


-03 


5400 


N&S 


720 


0.88 


23.61 


213 


06D3CC 


14:17:31.604 


+52:54:44.79 


2006 


-02 


-02 


5400 


N&S 


720 


0.58 


23.55 


99 


06D3cl 


14:22:15.106 


+52:16:41.53 


2006 


-03 


-06 


5400 


N&S 


720 


0.61 


23.75 


264 


06D3cn 


14:19:25.826 


+52:38:27.78 


2006 


-03 


-06 


1800 


C 


680 


0.40 


21.26 


372 


06D3du 


14:21:18.563 


+53:03:27.10 


2006 


-05 


-24 


5400 


N&S 


720 


0.45 


24.27 


54 
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Table 6. continued. 



SN Name 


RA (2000) 


Dec (2000) 


UT Date 


Exposure (s) 


Mode 


•^cent 


seeing (") Mag (i') 


%inc 


06D3dz 


14:18:14.731 


+52:56:23.11 


2006-05-02 


2400 


C 


680 


0.98 22.19 


66 
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Table 7. SNLS identifications and derived properties for SNe candidates observed from Nov. 2004 to May 2006. 



SN z 








CI & Type" 


Template'^ 


Z from" 


04Dlpu 


0.639 


0.001 


3 - 


SN la* 




on 


04D2mh 


0.59 


0.01 


3 - 


SN la* 


1996X +2 


SN 


04D2mj 


0.513 


0.001 


5 - 


SN la 


1992A-1 


on, H/3, OIII 


OSDlaz 


0.842 


0.001 


5- 


SN la 


1996X +2 


possible on 


OSDlby 


0.298 


0.001 


5- 


SNIa 


2002bo-l 


Hot, Hj8, OIII, SII 


OSDlbz 






2- 


SN 






OSDlcc 


0.564 


0.001 


5 - 


SN la 


1996X +7 


HyS, on, OIII 


OSDlee 


0.558 


0.001 


4- 


SN la 


1992 A +6 


H&K 


05Dlej 


0.312 
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Table 7. continued. 
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SN z (tI CI & Type'' Template" zfrom" 

a — redshifts identified with host galaxy features are accurate to » +0.001; redshifts identified 
with template matches are accurate to « ±0.01 

b — final object identification and confidence index 
c — best low-z spectral fit (name and days past max.) 
d — the host features used to estimate z 
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Appendix A: Results - Spectroscopic analysis of 
SNLS SNe la 
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Table A.3. Equivalent width and ejection velocity results for the published high redshift SNe la observed by the SNLS at the Gemini telescopes. 
This set includes all of the SNLS objects (reduced by the author) that were confirmed as Type la SNe (a confidence index of 3 or higher) and 
were subject to less than 65% contamination from their host galaxies. The blank spaces ('■ ■ ■') indicate where a measurement had to be removed 
due to inadequate wavelength coverage. The uncertainties in the EW results include measurement errors, errors from any possible variance in the 
pseudo-continuum, and the uncertainty from the host galaxy contamination correction. 
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Name Day^'^ EWjCalTF EWjSill}'^ EW{Mglli'" \'^' 7 1000.0 

El 

a — Rest- frame day, relative to B' max 

b — Equivalent Width, A 

c — Ejection velocity, km/s 

* Icr errors listed in parentheses 
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Appendix B: SNLS Observed Spectra 

Sprectra from each Supernova Legacy Survey SN candidate ob- 
served at the Gemini telescopes from December 2004 through 
May 2006. The observed properties and observation setup for 
each object is summarized in Tableland the final object iden- 
tifications and redshift estimates are displayed in Table|7] These 
spectra have been corrected for sky absorption, the narrow host 
galaxy features have been removed, and the GMOS chip gaps 
have been interpolated over. These images show the observed 
spectrum, re-binned to 7.5 A resolution, in dark blue with a 
smoothed version of the spectrum overplot in light blue. The 
smoothed spectrum is for illustrative purposes only and was gen- 
erated witha Savitsky-Golay filter computed with a 3"^ degree 
polynomial fit over 40 pixels of the observed spectrum. Note 
that for some of these spectra a large amount of host galaxy con- 
tamination is clear. 
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Rest Wavelength [A] Rest Wavelength [A] Rest Wavelength [A] 

3000 3500 4000 4500 5000 5500 6000 3000 3500 4000 4500 5000 5500 6000 3000 3500 4000 4500 5000 5500 6000 6500 




6000 7000 6000 9000 6000 7000 6000 9000 6000 7000 BOOO 9000 

Observed Wavelength [A] Observed Wavelength [A] Observed Wavelength [A] 




6000 7000 6000 9000 5000 6000 7000 8000 6000 7000 3000 9000 

Observed Wavelength [A] Observed Wavelength [A] Observed Wavelength [A] 



Rest Wavelength [A] Rest Wavelength [A] Rest Wavelength [A] 

2500 3000 3500 4000 4^00 5000 5500 4000 5000 6000 7000 BOOO 3500 3000 3500 4000 4500 5000 5500 




6000 7000 BOOO 9000 5000 6000 7000 8000 5000 6000 7000 BOOO 

Observed Wavelength [A] Observed Wavelength [A] Observed Wavelength [A] 




6000 7000 BOOO 9000 5000 6000 7000 BOOO 6000 7000 8000 9000 

Obb-LTVcd WavL-lciigth [A] Obsc-rvcd WavL-lciigth [A] Observed Wavclungth [A] 



Fig. B.l. Confirmed SNe la. The observed, re-binned spectrum is shown in dark blue and the smoothed spectrum is overplot in light 
blue for illustrative purposes. 
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Rest Wavelength [A] Rest Wavelength [A] Rest Wavelength [A] 

2500 3000 3500 4000 4500 5000 3000 3500 4000 4500 5000 5500 6000 3000 4000 5000 6000 7000 




eooo 7000 Booo 9000 eooo 7000 booo 9000 6000 7000 eooo 9000 

Observed Wavelength [A] Observed Wavelength [A] Observed Wavelength [A] 



Rest Wavelength [A] Rest Wavelength [A] Rest Wavelength [A] 

3000 3500 4000 4500 5000 5500 6000 3000 4000 5000 0000 7000 3500 3000 3500 4000 4500 5000 5500 




6000 7000 8000 9000 5000 6000 7000 6000 5000 6000 7000 BOOO 

Observed Wavelength [A] Observed Wavelength [A] Observed Wavelength [A] 




Fig. B.2. Confirmed SNe la. The observed, re-binned spectrum is shown in dark blue and the smoothed spectrum is overplot in light 
blue for illustrative purposes. 
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Fig. B.3. Confirmed SNe la. The observed, re-binned spectrum is shown in dark blue and the smoothed spectrum is overplot in light 
blue for illustrative purposes. 
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Fig. B.4. Unidentified SN candidates. The observed, re-binned spectrum is shown in dark blue and the smoothed spectrum is overplot 
in light blue for illustrative purposes. 



